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Abstract 

The global analysis of S ^ PP, PV charmless decays with QCD factorization (QCDF) is 
presented. The predictions of QCDF are in good agreement with experiments. The best fitted 
CKM angel 7 is around 79°. The predicted branching ratios oi B ^ tt^tt^, ujK~^ , wtt"*", t:^K*^ etc. 
are also in good agreement with new data of BaBar and Belle. 
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I. INTRODUCTION 



The charmless two-body B decays are very important for extracting CKM angles a, 7 
and for testing QCD. Up to now, BaBar and Belle have accumulated large set of data. It 
is highly interesting to analysis these data by using different theories and compare the data 
with theoretical predictions. QCD factorization [l| has been used to analysis B PP, PV 
charmless decay data (here P denotes pseudoscalar meson, V vector meson). With the 
data at that time, the theoretical results prefer a larger CKM angle 7. The QCD factorization 
(QCDF) predictions for some B PV channels are only marginally consistent with the 
experimental data. Notice that the QCDF predictions contain large numerical uncertainties 
due to the CKM matrix elements, form factors, and annihilation parameters. Furthermore, 
the uncertainties of various decay channels are strongly correlated to each other. So we are 
stimulated to do a global analysis to check the consistency between the QCDF predictions 
and the experimental data. Beneke et al. ^ have done a gloabl analysis for B —>■ PP, 
including tttc, tcK modes with QCDF approach. The results show a satisfactory agreement 
between QCDF predictions and data for B PP branching fractions. But their fitted 
CKM angle 7 ~ 90° which is not consistent with the standard CKM fit 5]: 37° < 7 < 80°. 
Now BaBar and Belle have a lot of data on B ^ PV. It is necessary to do the global 
fit to i? ^ PP and PV at the same time. We have done it and found that QCDF can fit 
all data on B ^ PP, PV. The fitted angle 7 ~ 79° which is consistent with the standard 
CKM global fit. We also have new predictions on i? — ^ 

My talk is organized as follows: in SectionlHl I give general remarks on QCD factorization. 
Section UTTl is devoted to the global analysis on i? ^ PP, PV decays. In Section HV| I give 
the summary and conclusion. 

II. GENERAL REMARKS ON QCD FACTORIZATION 

The low energy effective Hamiltonian for = 1 is 
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The decay amplitude for B M1M2 is 

AiB^AhM^) = %Y.Y.C^mM,M,m^^)\B), (2) 

i q 

The Wilson coefficients Cj(yu) include short distance effects above the scale 0{mb) and 
are perturbatively calculable. The long distance effects (contributions from the scale below 
(9(mb)) are included in the hadronic matrix elements {MiM2\Qi{fi)\B) . 

For naive factorization (BSW), 

(MiM2|Q.(/^)|B) ~ (Mi|Ji|0)(M2|J2|5) 

y \ (3) 

decay constant form factor 
There are no renormalization scheme and scale dependences. So the corresponding de- 
pendences of connot be canceled and the "non-factorizable" contributions cannot be 
accounted for. 

In 1999, Beneke et al. jl| proposed a new factorization scheme based on QCD. That is 
so-called QCD factorization (QCDF). In this scheme, in the heavy quark limit, 

+ E j\difdxj\Tll{i,x,y)^B{i)^MA^)^MM+0{KQCD/m,) (4) 

where T/'^^ denote hard scattering kernels which are calculable order by order in perturbation 
theory and include short distance effects. The long distance effects are included in decay 
constants, form factors and ling-cone distribution amplitudes. Mi is a light meson or a 
charmonium state. M2 (containts the spectator quark in B meson) is any light or heavy 
meson. If M2 is heavy (for example, D meson), the second line in the Eq.(0]) is 1/m;, 
suppressed. 

General observation: 

• At the zeroth order of a^, it reduces to naive factorization 

• At the higer order, the corrections can be computed systematically 

• The renormalization scheme and scale dependence of {Qi) is restored. In the heavy 
quark limit, the "non-factorization" contributions is calculable perturbatively. It does 
not need to introduce Nf^ 
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FIG. 1: T^'^^ at 0{as). The upward quark lines denote the ejected meson Mi 




FIG. 2: annihilation contributions 

• In heavy quark limit, strong phases arise solely from vertex and penguin corrections, 
so are at 0{as). That means that strong phases are a g- suppressed. 

• Numerically, KqcD/i^h ~ ^{.(^s)) so power corrections contribution to strong phases 
are also important. We do not have a systematic way to estimate power {h^Qco/^b) 
corrections. Thus the calculation of strong phase is not reliable. So CP asymmetry 
calculation is unreliable. 

• ly-exchange & VT-annihilation diagrams are l/m^ suppressed 

• No long distance interactions between Mi and {BM2) 

• includes: 

a° : tree diagram 

as '. non-factorizable gluon exchange 
vertex corrections (Fig|Tl^a)-(d)) 
penguin corrections (Figme)-(f)) 
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• T includes: 

as : hard spectator scattering (Figmg)-(h)) 

Superficially, KqcD/f^b ~ 1/15 is a small number, But in some cases, such power suppression 
fails numerically. For example 



is always multiplied by a formally power suppressed but chirally enhanced factor = 
2fip/mb {fip = mp/(mi + 7712), rrii are current quark masses). In the heavy quark limit, rrib 
^ 00, ^ 0. But for nif, ~ 5GeV, ~ C^(l)) no I /nib suppression. So we should consider 
the chirally enhanced power corrections. For example, in B ^ tcK, penguin is important. 
Dominant contribution to the amplitude ~ 04 + agr^, (04 ~ dQr^). So chirally enhanced 
power correction term a^r^ cannot be neglected. 
Possible sources of power corrections: 

• high twist wave functions; 

• quark transverse momentum k±; 

• annihilation topology diagrams. 

When we include chirally enhanced power corrections, the infrared (including soft and 
collinear) divergences for vertex corrections cancel only if $o-(a^) is symmetric, i.e. ^a{x) = 
<l>o-(l ~ ($0-, $p are twist-3 wave functions). But even $o-(x) = $o-(l — x), there is still 
logarithmic divergence for hard spectator scattering and annihilation topology. It violates 
factorization. 



Other power corrections maybe unimportant based on Renormalon Calculus estimation 

III. GLOBAL ANALYSIS OF B ^ PP, PV 

The decay amplitudes of i? ^ PP, PV can be written as 




(5) 




A/,~0.5GeV, < ph,a < 1, 0°<(f)H,A < 360° 



(6) 
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where hi = bi{XA,<ys) = o^sfiXA) is related to the contributions of annihilation topology. 

We use CKMFlTTER package [3] developed for global analysis of i? ^ PP and enlarged it 
to include B PV. The RFiT scheme is implemented for statistical treatment, assuming 
the experimental errors to be purely Gaussian, while quantites that cannot be computed 
precisely are bound to remain within their predefined allowed ranges. 

The input parameters are as follows: 

1. two-parton light-cone distribution amplitudes. 



for pseudoscalar meson 

{P{k)\q^{z^)qp{z,)\0) 
ifp 



['dxe'^''-^^+''-'^nH<!>p{x) -^^p^,UPp{x)-a,,k^z''^^]}^ (8) 

JO L 5 J J /3o 
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where z = zi — Z2. 

for longitudinally polarized vector meson 

for vector meson, the contributions of the twist-3 LCDAs are doubly suppressed, 
so can be safely disregarded. 

In order to reduce the number of input parameters, we use the asymptotic light- 
cone distribution amplitudes for final light mesons instead of Gegenbauer polyno- 
mials. 

$p(x) = Qxx, $p(x) = 6xx, = 1, ^\{x) = 6xx, (10) 

for B meson, its wave function appears only in the contributions of hard spectator 
scattering jj| 

niB 



A 



B 



Xb = (350±150) MeV, (11) 
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2. chirally enhanced factor & annihilation- related parameters 

2777,2 

r"~r:^~rf = ^, mJmfe) = 4.2 GeV, m J2 GeV) = (110±25) MeV, 

^ ^ ^ nismb 

XZ = ln^(l + pZe'^"'^^)^ = ln^(l + pZ^'''"-) (12) 

(note: X^, for B ^ PP are different for 5 ^ Pl^) 

3. Decay constants and form factors 

= 131 MeV, Jk = 160 MeV, Jk* = 214 MeV, /p = 210 MeV, 
/.^ = 195 MeV, = 233 MeV, /b = 180 MeV, (p = 39.3°, 
/, = 1.07^, /, = (1.34±0.06)^, /^^ = /,cos0, = -/.sin0, 

F^^p'^(O) = 0.28±0.05, A^"^(0) = 0.30±0.05, R^k=¥^b^ = 0.9±0.1, (13) 



■/if 



f,^\v,) = (0|s75s1r/) = (/,^ - ^8 = 0- arctan(V2/,//,), 

= ks), tefM = |, ^0 = 0- arctan(v^/,//,), 

4. CKM parameters 

We use Wolfenstein parameterization for the CKM matrix and take 

A = 0.835, A = 0.22, IK^I = 3.49±0.24exp±0.55theo, (14) 

and use the measured branching fractions listed in Tablefljas input. We did not use the 
branching ratios in TablelTTlas input. The reason is that the errors of the data are large 
and the glue content of r/^'^-meson cannot be treated neatly. Detailed discussion can 
be found in our published paper j^. For experimental constraints of CP asymmetries, 
the corresponding QCDF predictions are not reliable. Thus we do not use measured 
CP asymmetries as input in our global fit. 

For B TTTT, TxK, the fit is sensitive to \Vub\, l{p,v)^ pB^-K ^ F^^^ , Xa, /b/^b, 
Including seven B^PV decay channels, only Aq^'', X^ are newly involved sensitive pa- 
rameters. So including B —>■ PV we shall have more stringent test of QCDF. Our best fit 
of 7 ~ 79° (see FigureOI) which is consistent with recent fit results 37° < 7 < 80° Q]. The 
best fitted results are listed in Table UTTl 

The confidence levels of the fitted branching ratios are presented in FigurelU From 
FigurelUwe can see that 
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7 [degrees] 



FIG. 3: left: the contours indicate >90% C.L.& >5% C.L. right: C.L. of angle 7. 

• For B — i> 7r°7r°, the best fit is around lxlO~^ while the BaBar and Belle averaged 
measurement is (1.90±0.49)xl0~^. 

• For —>■ ujK~^, the best fit is 6.25x10^^ while the Belle measurement is 
(9.2^2:^±1.0)xl0^^ and the CLEO measurement is < 8xl0~^ 

• For B^ —>■ uJTT^, the best fit is 6.66x10"^ while the BaBar measurement is 
(6.6ti:8±0.7)xl0~^ and the Belle measurement is < 8.2x10^^ 

• For B^ 71^ K*^, the best fit ~ lOx 10^^ while the BaBar preliminary measurement 
reported at LP03 is (10.3±1.2t^:?)xl0-^ 

On the whole, our global fit is successful. 



IV. SUMMARY & CONCLUSION 

• QCD factorization is a promising method for charmless two-body B decays, which are 
cruical for the determination of the unitarity triangle. 

• We enlarged the CKMFitter package to include B — * PV charmless decay channels 
and did a global analysis. It is shown that the QCDF predictions are basically in good 
agreement with the experiments. 

• We obtain 7 ~ 79°, consistent with the CKM global fit. 
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FIG. 4: The confidence levels of the fitted branching fractions (in the unit of 10 ^). 

• For B 7r°7r°, the best fit is around lxlO~^ while the BaBar and Belle averaged 
measurement is (1.90±0.49) xlO~^. 

• For Oi)K~^, the best fit is 6.25 xlO~^ while the Belle measurement is 
(9.2l|^±1.0)xl0-^ and the CLEO measurement is < SxlQ-*'. 

• For S+ — > a;7r+, the best fit is 6.66x10"^ while the BaBar measurement is 
(6.6l?:8±0.7)xl0-^, and the Belle measurement is < 8.2 xlO"^. 
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• For TT'^K*'^, the best fit ~ 10x10 ^ while the BaBar prehminary measurement 

reported at LP03 is (10.3±1.2l^:?) x 10"^ 

The results are already published in 0|. I thank Drs. Junfeng Sun, Deshan Yang and 
Guohuai Zhu for their collaboration. 
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TABLE I: Experimental data used in the global fit 



HrxlO^ 


CLEO 


BaBar 


Belle 


weighted 
average 




4.3+|i±0.5 

1.41: 


4.7±0.6±0.2 


5.4±1.2±0.5 


4.77±0.54 


B+ ^ tt+ttO 


— z.yj 


5.5+iS±0.6 

— u.y 


7.4+^-^±0.9 


5.78±0.95 


K+TT- 


17.2+^i±1.2 


17.9±0.9±0.7 


22.5±1.9±1.8 


18.5±1.0 


B+ ^ i^+7rO 


— Z. i — l.O 


12.8+H±1.0 

— 1 . ± 


13.0+2i±l-3 

— Z.4 


12.7±1.2 


S+ ^ K07r+ 


18.2+!-^±1.6 


17.5+H±1.3 

— 1. / 


19.4+H±1.6 

— O.U 


18.1±1.7 




— O.i — o.o 


10.4±1.5±0.8 


8.0+|?±1.6 


10.2±1.5 


5+ -> 7?7r+ 


< 5.7 


< 5.2 


5.3+?-? (< 8.2) 


< 5.2 




27.6+?1±4.2 

— / .4 


28.9±5.4±4.3 


on O+6.0+2.8 

— O.O — O.i 


25.4±4.3 




10.4+H±2.1 


< 39 


8.0+2n±0.7 


8.6ib2.0 


so ^ ^+^- 


16.0+Ii±2.8 

— 0.4 




-,1 9+5.9+1.9 

J- -L c K 1 o 
— O.O — i.o 


13.1±4.7 


B+ 4>K+ 


5.5l?i±0.6 


9.2±1.0±0.8 


10.7±1.0t?| 


8.9±1.0 




5.4t|^±0.7 


8.7tJ:^±0.9 




8.6±1.3 


B+ rip+ 


< 10 




< 6.2 


< 6.2 




< 21 


5.9li:^±o.9 




5.9±1.9 


TABLE IL data not used in the global fit 


Br xlO^ 


CLEO 


BaBar 


Belle 


weighted 
average 


B+ r]K+ 


< 6.9 


< 6.4 


5.2t\-l (< 7.7) 


< 6.4 


B+ Tr+K*° 


< 16 


15.5±3.4±1.8 


16.2l|^±2.4 


15.8±3.0 


B° Tr-K*+ 


16l^±2 




26.0±8.3±3.5 


19.0±4.9 


B+ riK*+ 


26.4:tli±3.3 


22. 11:^^2^ ±3.3 


26.5^^:^±3.0 


25.4±5.3 


BO riK*^ 


13.8l|J±1.6 


19.8l|^±1.7 


16.51^2^1-2 


16.4±3.0 


B+ — ^-/v ^ 


< 8 


< I 


9.2;^:^±i.o 




B+ uTr+ 


11.3l|^±1.5 




< 8.2 
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TABLE III: Fitl and Fit2 mean the best fit value with and without the contribution of the chirally 
enhanced hard spectator and annihilation topology, respectively. 



modes 


Exp. 


Fitl 


Fit2 


5° 


— > tt'^tt 


4.77±0.54 


4.82 


5.68 


B+ 




5.78±0.95 


5.35 


3.25 


5° 




18.5±1.0 


19.0 


18.8 


B+ 


-> K+TT^ 


12.7±1.2 


11.4 


12.6 


B+ 


K^Tr+ 


18.1±1.7 


20.1 


20.2 


5° 




10.2±1.5 


8.2 


7.3 


5+ 


— ^ 777r"'" 


< 5.2 


2.8 


1.8 


5° 


— y TT'^ 


25.4±4.3 


26.7 


29.5 


B+ 




8.6±2.0 


8.9 


8.5 


5° 


^ K+p- 


13.1±4.7 


12.1 


5.1 


B+ 


(f)K+ 


8.9±1.0 


8.9 


7.1 


5° 


^ (t>K^ 


8.6±1.3 


8.4 


6.7 


B+ 


rip+ 


< 6.2 


4.6 


3.8 


5° 




5.9±1.9 


6.3 


1.2 
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